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Abstract 

We study the reflection of a Hermite-Gaussian beam at an inter- 
face between two dielectric media. We show that unlike Laguerre- 
Gaussian beams, Hermite-Gaussian beams undergo no significant dis- 
tortion upon reflection. We report Goos-Hanchen shift for all the spots 
of a higher order Hermite-Gaussian beam near the critical angle. The 
shift is shown to be insignificant away from the critical angle. The cal- 
culations are carried out neglecting the longitudinal component along 
the direction of propagation for a spatially finite, s-polarized, full 3-d 
vector beam. We briefly discuss the difficulties associated with the 
paraxial approximation pertaining to a vector gaussian beam. 

1 Introduction 

The Goos-Hanchen (GH) shift, first discovered experimentally in 1947, 
refers to the longitudinal displacement of a beam from its expected 
geometric position under total internal reflection Artmann was 
the first to provide a theoretical explanation of the effect using the 
stationary phase approximation. A somewhat different calculation of 
the Goos-Hanchen shift using the conservation of energy was put forth 
by Renard Quantitatively the GH shift can be interpreted as the 
shift of the peak intensity spot of the beam [2] or the weighted mean 
position of the beam [4]. A thorough theoretical analysis of the effect 
was carried out by Lotsch [5] in 1968. The effect, although discovered 
quite a while ago, still remains a topic of interest in current optics 
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literature. The current interest in GH shift is further fueled by sev- 
eral means for enhancing the effect and its use for sensing [6]. These 
include structures involving absorbing media [71 [8], metals [9], left- 
handed metamaterials [inilll]- Very recent studies focus on systems 
supporting surface plasmon resonances [U [121 [E] 8.nd frustrated total 
internal reflection [TJ]. From a theoretical angle, most of the studies 
on GH shift used to suffer from two major drawbacks. First, an over- 
whelming majority deal with two dimensional fundamental beams (the 
third dimension is supressed) [151 [HI [IT] , and second, they overlook 
the vector characteristics of the beam. Recall that the paraxial elec- 
tromagnetic wave is not an exact solution of the Maxwell's equations, 
but holds only to the first order in ^ {k being the magnitude of the 
wavevector). Since the pioneering work by Agrawal and Pattanayak 
[18j . the difficulties associated with the paraxial approximation are 
now clear. The vector nature of the beams adds to the complexity 
of the problem (leading to a longitudinal field component). Only re- 
cently vector Gaussian beams have been dealt with in the context of 
propagation in a uniform medium [19] . A systematic treatment of GH 
shift for a fundamental vector Gaussian beam was given by Aiello and 
Woerdman [20j. It was shown by the same authors in a later paper 
that the suppression of the correction due to the longitudinal com- 
ponent recovers Artman formula for the shift There has been a 
great deal of interest in the Laguerre-Gaussian beams and their lon- 
gitudinal and lateral shifts in view of the angular momentum carried 
by the beam [211 [22l [23l [HI [SSI [26] . The longitudinal and transverse 
shifts for the Laguerre-Gaussian beams were calculated by Bliokh et al 
[24j and verified experimentally by Merano et al [25]. Another recent 
paper presents both theoretical and experimental results on the shift 
of Laguerre-Gaussian beams and its distortion under refiection and 
transmission [26]. However, no such studies exist for simpler Hermite- 
Gaussian beams. In this paper we concentrate on a general vector 
Hermite-Gaussian beam with s-polarization and look at the reflection 
of the same from a denser to lighter medium interface. We show how 
the longitudinal component arises from a simple description of the 
vector field by means of scalar and vector potentials. Neglecting the 
longitudinal component under the strong paraxial approximation, we 
use a angular spectrum decomposition to calculate the refiected field 
for each spatial harmonic. Summing up the harmonics we end up with 
the reflected beam profile. 

The structure of the paper is as follows. In section 2, we provide 



2 



a simple mathematical formulation of the problem and show that a 
longitudinal component is inevitable in the paraxial formulation. In 
section 3, we first validate our code by matching our numerical results 
with those predicted by Artman formula. Next, we present the con- 
tour plots for fundamental and higher order Hermite-Gaussian beam 
profiles. We also comment on the possibility of splitting of the re- 
flected beam, which has been observed for tightly focussed beams. 
Finally, we summarize the main results in the conclusion. 



2 Mathematical formulation 

Consider the electromagnetic wave of frequency to described by the 
electric E and magnetic B fields in a medium with dielectric constant 
e and permeability constant ^ = 1. Let the scalar and vector potentials 
be (j) and A, so that 

B = V X A = H, (1) 
1 8 A 

B = -V(^--— = -V^ + ikoA, (2) 
c ot 

where /cq = 7 is the vacuum wave vector. Under the Lorentz condition, 
the vector A satisfies the following equations: 

V • A - ikoect) = 0, (3) 

V^A + k^eA = 0. (4) 

One way to tackle the problem is to ignore the variation of all quanti- 
ties in one direction, making it effectively a two dimensional problem. 
Polarization is no longer an issue in the two dimensional scenario and 
the standard scalar approximation techniques can be used. A more 
general set of solutions is obtained by assuming that the vector A, 
rather than E is polarized in a particular direction. 

We seek the solution of A in the form of a wave propagating along 
z with transverse polarization (say along x) . Such a solution takes the 
form: 

A = [a(x,y,2)e*'''^ + c.c]x, (5) 

where k = k^y/e and a{x, y, z) is the slowly varying (with respect to 
z) amplitude. Using equations ^ and ([5]), Vcj) can be evaluated to be 

VcP = iiko€)-^V exp (ikz)^ , (6) 
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which subsequently leads to the expression of E as follows. 

E = iko (ax + — — — z | exp (ikz). (7) 
\ k ox J 

It is clear from equation ([7]) that the paraxial approximation inevitably 
leads to a longitudinal component. Substituting expression ([5]) in 
equation (j4]) and making use of the slowly varying approximation one 
can write down the paraxial equation for a light beam: 



(8) 



d'^a d'^a da 
oy-^ oz 

The solution of ([S]) can be expressed as: 

aix,y,z) = J a{q.j;,qy)ex.p[i{q^x + qyy)]exp (^-i ^ ^"^ ^ z^dq^dqy. 

(9) 

For localized profiles of the form a{qx,qy), one has solutions in the 
form of a beam. For example, for 

a{q.,qy) = vrT^o' exp ^-Wo' ^Sl±3l^^ , (10) 

one recovers the fundamental gaussian beam: 

(. 



Wo ( —r^ \ ( ikr^ 



(11) 

The expressions for the beam waist W{z), radius of curvature R{z) and 
Gouy phase (p{z) in the above expression are standard. Following a 
similar route, the expression for the general Hermite-Gaussian beam 
can be calculated. In order to calculate the reflected beam profiles 
for higher order Hermite-Gaussian modes, we have taken the incident 
beam amplitude as follows: 

am,n = irfeexp(^^)exp(^)exp(-i(^(z)„,„) 
xexp{ikz)K^(^)HJ^), 



w(z) J V'^(-^) 

where and R{z) have the same meaning as before, but the Guoy 

phase takes on the definition: 



f m,n{^) = {m + n + 1) tan ^ (^-^^ , (13) 



where Zr = 

In the context of the problem of calculating the reflection coef- 
ficient, we now specify the geometry. The geometry is akin to that 
considered by Aiello and Woerdman [1]. We consider an interface (at 
z = 0) between vacuum and a dielectric medium with dielectric con- 
stant e. Let the beam be incident at an angle to z and with vector 
potential polarized along y. Let also z' axis coincide with the direc- 
tion of propagation of the beam. The primed and unprimed axes are 
oriented at 6 to each other. The relation between the unprimed and 
primed co-ordinates are given by 

x' = xcosO — zsm.9, y' = y, z' = zcosO + xs'mO. (14) 

For a fundamental gaussian beam of the form given by (jlOp . E reads 
as follows. 



E = y-^z' a(,i 



X exp 



' ,2,2 N 

QxX +qyy ^ 



d'q'. (15) 



This expression is consistent with the result of Chen et al [19j, that 
the choice of a zero transverse field component at the input boundary 
plane ensures the vanishing of the component throughout the prop- 
agation in the same medium. The consistency of the magnitude of 
the total wave vector of every spatial component under the paraxial 
approximation \J(i^^\-qy «k can easily be verified. 

e + <iy' + k\l-'^^^f^k\ (16) 

Note also that under very weak transverse variation of the field when 

the inequality -^{« 1) holds, ([T5|) leads to a simple s-polarized beam 
all throughout [26j. Henceforth, we refer to this approximation as the 
strong paraxial approximation. 

Moving to the unprimed coordinate system one can write down a 
typical plane wave component of E as: 

(cj^x + ei^y + ei^z)a{q^,qy) exp [i{qxX + Qyy + q^z)], (17) 
where e^^ = -^sinO, e^^ = 1, Cj^ = -^cos9 and qx, qy and qz are 
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given by: 

Qx = Qx cos 9 + ksinO I 1 



/ 2 I ,2- 
Qx +% 

2P 



Qy = dy, (18) 



qz = q'^sinO — k cos 9 ( 1 



/ 2 ,2' 
ix q.y 

2A;2 



It is again easy to check the consistency of the total wave vector for 
a spatial component under the paraxial approximation using above 
equations 

Qx^ + qy^ + Qz^ « k"^- (19) 

In what follows, we neglect the terms proportional to | in the 
amplitude, while retaining the full dependence in the phase factors. 
Thus the beam reduces to an s-polarized beam with only non- vanishing 
component along y. Then, one can easily calculate the reflected com- 
ponent by multiplying the amplitude of the spectral component with 
the reflection coefficient which is given by 

1 _ 9£a 



where Qz^ = \/k'o — {qx"^ + Qy^)- In order to validate the above approx- 
imation, one has to have a broad beam (close to plane wave profile), 

so that " << 1. The broader the beam, the smaller is the GH 
effect. Note that plane waves do not exhibit any GH shift. Thus one 
has to exercise care in order to choose the parameter domain where 
the shift, albeit small, is still observable, if one wants to describe 
the phenomenon with a model of an incident gaussian beam with a 
given polarization. For a tightly focused beam, the paraxial approx- 
imation forces one to include the longitudinal component. With the 
longitudinal component present, the calculation of the reflection and 
transmission coefficients poses a formidable problem. 



3 Numerical results and plots 

In this section we follow the procedure detailed in the previous section 

to calculate the beam profile on the interface before and after the re- 
fiection. To be specific, we take a beam described by a scalar function 
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with s-polarization. We find the 2-d Fourier components incorporat- 
ing the obhque incidence. Each Fourier component is multiphed by 
the corresponding reflection coefficient. Then an inverse Fourier trans- 
form leads to the reflected beam profile at z = 0. We have considered 
the cases of fundamental and higher order Gaussian beams. We show 
that identity of the beam is retained in the case of Gaussian-Hermite 
beams. In all our calculations we have taken ^/€ = 1.7321 (the medium 
of incidence), eo = 1.0, kWo = 40, unless specified otherwise. The in- 
cident beam is assumed to be focused at z = (at the interface). 

While choosing the relevant parameters, special attention was given 
to two aspects. First, the sampling frequency was taken to be more 
than the Nyquist rate |27j . Second, the parameter domain was chosen 
such that the strong paraxial approximation was valid. For an angle 
of incidence of 9 = 40° and kWo = 40, for the y direction, we take the 
sampling interval to be Ay/Wo = 0.1020. Sampling is done inside a 
square area of side of 15 units. In the inverse Fourier space, the sam- 
pling interval is Aqy = 0.0602. Let Qy represent the value at which 
the amplitude falls off, say, to 10~^ times the peak value. The value 
of Qy gives an estimate of the angular dispersion of the incident beam. 
For a fundamental Gaussian beam with kWo = 40, we have qy ~ 0.17. 
Thus the chosen parameters alleviate both the above concerns. 

3| , , n • • 1 
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Figure 1: Normalized GH shift as a function of angle of incidence 9. Solid 
line: numerically computed GH shift; dashed line: shift predicted from the 
Artmann formula 

We first use the numerical method to predict the GH shift for a two 
dimensional Gaussian beam. Figured] shows the dependence of the GH 
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shift (normalized to the beam waist Wq) on the angle of incidence 9 
for a 2-d fundamental Gaussian beam. The GH shift values computed 
using the procedure described in the previous section (solid line) are 
contrasted with those predicted by the Artmann formula [15] (dashed 
line) . It is clear from figure [J that for angles larger than the critical 
angle and away from it, the GH shift computed numerically falls off 
and compares well with that predicted by the Artmann formula. Near 
the critical angle, it is well known that the Artmann formula has a 
singularity [28]. Thus the stationary phase approximation, which is 
used to derive the Artmann formula, is not applicable near the critical 
angle. Our numerical results are consistent with the fact that the GH 
shift must show a smooth, continuous and finite variation near the 
critical angle [28j . 

Having validated our numerical procedure in the two dimensional 
case, we now use it to analyze three dimensional Gaussian beams. 
The profiles on the interface for the incident and reflected fundamen- 
tal Gaussian beam are shown in figure [2j The left (right) pane in each 
horizontal line shows the contour of the incident (reflected) beam at 
the interface. Horizontal (vertical) axis is the normalized x{y) di- 
rection. It is clear from figure [2] that for small and large angles of 
incidence, there is practically no shift of the fundamental Gaussian 
beam. Close to the critical angle 9cr (= arcsin(Y^eo/e)), the shift is 
discernible. Slightly above the critical angle, the shift achieves the 
maximum value. 

Figure [3] shows the results for TEMu beam incident at 40° (ap- 
proximately where the GH shift is largest). It is clear from Figure 
[3] that even after reflection the beam identity of the TEM beam is 
intact, though there is a slight difference in the heights of the pair 
of lobes. The shift becomes negligible for small and large values of 
the angle of incidence. Figure [3] also shows the incident and reflected 
profiles for TEM21 (middle) and TEM22 (bottom) cases, respectively. 
Even for For Hermite-Gaussian beams the GH shift is small for small 
and large angles of incidence, in this case peaking at around 40°. The 
GH shift for near normal incidence is almost negligible. 

We now focus our attention to beam distortion under reflection 
near the critical angle. For loosely focussed beams, it was shown that 
the distortion is negligible for Hermite-Gaussian beams. This is not 
the case for Laguerre- Gaussian beams because of the angular momen- 
tum carried by the beam. The inherent orbital angular momentum 
|22j of the beam and the lateral Feodorov-Imbert shift are responsi- 
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Figure 2: Incident (left) and reflected (right) profiles for fundamental Gaus- 
sian beams. The angles of incidence are top: 10°, middle: 40°, bottom: 
60° 
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Figure 3: Incident (left) and reflected (right) beam profiles for TEMn (top), 
TEM21 (middle), TEM22 (bottom), incident at 40°. 
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ble for the distortion of the beam profile upon reflection [26]. Thus, 
loosely focused higher order Hermite-Gaussian beams can retain their 
identity upon reflection, while the Laguerre-Gaussian beams can be 
distorted beyond recognition. 
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Figure 4: The reflected beam profile for 2-d (left) and 3-d (right) Gaussian 
beam incident at 25°. 

Until now, we have been looking at broad (loosely focused) beams 
which fall well within the realm of validity of the paraxial approxima- 
tion. We now look at tightly focused beams which fall at the periphery 
of this regime. We qualitatively describe the effects of reflection on 
a tight Gaussian beam. Figure H] shows the incident and reflected 
profiles for a tightly focused {kWo = 15) fundamental beam for 2-d 
and 3-d simulations (left and right panels, respectively, in the figure). 
The refiected profile shows a splitting in the beam, with the newly 
developed lobe having very little intensity compared to the main lobe. 
This observation is in step with other investigations of a similar kind 
where beam splitting upon refiection has been reported [HI [29]. It 
has to be kept in mind that very tight beams, by default, fall in the 
regime of non-paraxial beams and that in this regime the longitudinal 
component cannot be summarily neglected. 



4 Conclusion 

In conclusion, we studied the refiection of Gaussian beams from a 
dielectric-vacuum interface and showed that for tight beams it is ad- 
visable to take into account the longitudinal component of the field un- 
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der paraxial approximation, while one can use standard scalar beams 
with some given polarization in the opposite case. We numerically 
studied the reflected beam profiles and demonstrated the GH shift 
for fundamental and higher order Gaussian beams. Calculations were 
carried out after neglecting the longitudinal component of the electric 
field, which is justified in the strong paraxial approximation regime. 
We also concluded that higher order Hermite-Gaussian beams retain 
their identity after reflection unlike Laguerre-Gaussian beams. For 
tight beams it is indicative that beam splitting occurs upon reflec- 
tion, but for more precise results the longitudinal component must 
be retained in this case. Further investigation is needed to tackle the 
problem of reflection and associated GH shift of vector beams with- 
out making the strong paraxial approximation. Work is underway on 
enhancing the GH shift of higher order Gaussian beams using various 
layered media with metal-dielectric structures and metamaterials. 
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